Dielectric behavior of oblate spheroidal particles: 
Application to erythrocytes suspensions 



J. P. Huang and K. W. Yu 

Department of Physics, The Chinese University of Hong Kong, Shatin, NT, Hong Kong 

Abstract 

We have investigated the effect of particle shape on the eletrorotation (ER) 
spectrum of living cells suspensions. In particular, we consider coated oblate 
spheroidal particles and present a theoretical study of ER based on the spec- 
tral representation theory. Analytic expressions for the characteristic fre- 
quency as well as the dispersion strength can be obtained, thus simplifying 
the fitting of experimental data on oblate spheroidal cells that abound in the 
literature. From the theoretical analysis, we find that the cell shape, coating 
as well as material parameters can change the ER spectrum. We demonstrate 
good agreement between our theoretical predictions and experimental data 
on human erthrocytes suspensions. 
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I. INTRODUCTION 



Dielectrophoresis and electrorotation (ER) offer a unique capability of monitoring the 
dielectric properties of dispersions of colloids and biological cells. Under the action of ex- 
ternal fields, these particles exhibit rich fluid-dynamic behaviors as well as various dielec- 
tric responses. It is of importance to investigate their frequency-dependent responses to 
ac electric fields, which yields valuable information on the structural (Maxwell- Wagner) 
polarization effects |T|^]. The polarization is characterized by a variety of characteristic 
frequency-dependent changes known as the dielectric dispersion. In particular, for the j3- 
dispersion (also called the Maxwell- Wagner dispersion, ranging from KHz to MHz), it may 
suffice to focus on the induced dipole moments of the particles. In this work, we will analyze 
the /3-dispersion based on the spectral representation theory. 

In the last two decades, various experimental tools have been developed to analyze the 
polarization of biological cells - dielectric spectroscopy || , dielectrophoresis ||] and ER 
techniques. Among these techniques, conventional dielectrophoresis and ER are usually 
applied to analyze the frequency dependence of translations and rotations of individual cells 
in an inhomogeneous and rotating electric field, respectively [|]||. Moreover, one is able to 
monitor the cell movements by using automated video analysis || as well as light scattering 
methods 0. As a matter of fact, the general cause of ER is the existence of a phase 
difference between the field-induced dipole moment and the external rotating field, resulting 
in a desired torque which causes the cells to rotate. 

In the dilute limit, the ER of individual cell can be predicted by ignoring the mutual 
interaction between the cells. However, the cells may aggregate under the influence of the 
external field. In this case the Brownian motion can be neglected, and the cell system 
becomes non-dilute even though it is initially dilute. As an initial model, we have recently 
studied the ER of two approaching spherical particles in the presence of a rotating electric 
field [0]. We showed that when the two particles approach and finally touch, the mutual 
polarization interaction between the particles leads to a change in the dipole moment of 
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individual particles and hence the ER spectrum, as compared to that of isolated particles, 
via the multiple image method ||. 

In this work, we consider further the effects of cell shape on the ER spectrum. In fact, 
a cell may readily deviate from a perfect spherical shape due to many reasons, e.g., under 



oblate spheroidal shape, such as human erythrocytes. 

Regarding the spectral representation approach [jlTjl , it is a rigorous mathematical for- 
malism of the effective dielectric constant of a two-phase composite material. Also, the 
approach was extended to deal with a three-phase material [0. The spectral approach of- 
fers the advantage of the separation of material parameters (namely the dielectric constant 
and conductivity) from the cell structure information, thus simplifying the study. From 
this approach, one can readily derive the dielectric dispersion spectrum, with the dispersion 
strength as well as the characteristic frequency being explicitly expressed in terms of the 
structure parameters and the materials parameters of the cell suspension ||13[ . In the present 
work, the analytic expression for the characteristic frequency is derived with the aid of the 
spectral representation approach, thus simplifying the fitting of experimental data. We show 
good agreement with experimental data on human erythrocytes JTjfl. 



Consider a suspension of spheroidal particles of complex dielectric constant e\ coated 
with a shell of e s dispersed in a host medium of €2, where 



with / being the frequency of the applied field Eo = Eqz. The depolarization factors of the 
spheroidal particles is described by a sum rule 



a hydrostatic pressure ||. In addition to a prolate spheroidal shape |Tt|, there exist cells of 



II. FORMALISM 
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where L z (l/3 < L z < 1) and L xy (— L x = L y ) are the depolarization factors along the z- and 
x- (or y-) axes of the oblate spheroidal particle, respectively. L z = L xy = 1/3 just indicates 
sphere, while L z — 1 disk. 

The phenomenon of ER is based on the interaction between a rotating electric field E 
and the induced dipole moment M. The dipole moment of the particle arises from the 
induced charges that accumulate at the interface of the particle. In circularly polarized 
fields, the axis of lowest polarizability, namely z— axis of the present particle, should be 



oriented perpendicular to the plane of field rotation |L5 . 

The angle between M and E is denoted by 6, where 9 = us x time and u = 2irf is 
angular velocity of the rotating electric field. The torque acting on the particle is given by 
the vector cross product between the electric field and the dipole moment, so that only the 
imaginary part of the dipole moment contributes to the ER response. In the steady state, 
the frequency-dependent rotation speed £l(f), which results from the balance between the 
torque and the viscous drag, is given by 

n(f) = -F(e 2 , E, r^/m^cos 2 6} + b y (sin 2 9)], 

= -F{e 2 , E, V )Im[b x ] = -F(e 2 , E, r))Im[b y ] = -F(e 2 , E, r))Im[b xy ], (3) 

where Im[- ■ •] indicates taking the imaginary parts of [•••], the angular brackets denote a 
time average, and F is a coefficient which is proportional to the square of magnitude of 
field E 2 , but being inversely proportional to the dynamic viscosity rj of the host medium. 
For spherical cells, F = e 2 E 2 /2rj. Regarding rj, spin friction expression suffices for spherical 
particle. However, for spheroidal particles or many particles interacting in a suspension, we 
must consider the more complicated suspension hydrodynamics. Since the angular velocity 
of the rotating field is much greater than the electrorotation angular velocity, i.e., ui ^> Q, 
the time averages are just equal to 1/2. For a single coated spheroidal particle, the dipole 



factor b z is given by [10.16 



£ _ 1 (gj ~ ^)[e s + L xy (ei - e s )} + (ei - e s )y[e s + L xy (e 2 - e s )] 

xy 3 (e a - ei)(e 2 - e a )yL xy {\ - L xy ) + [e 8 + (e t - e s )L xy \[e 2 + (e s - e 2 )L xy ] 
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where y is the volume ratio of the core to the whole coated spheroid. 

We are now in a position to represent b xy in the spectral representation. In what follows, 
we will show that from the spectral representation, we can obtain the analytic expression for 
the characteristic frequency at which the maximum ER velocity occurs. Let Z\ = £2(1 — 
and assume x = e s /e 2 , we obtain 

b xy = NP + j^— (5) 



where NP denotes the nonresonant part |L2] which vanishes in the limit of unshelled 
spheroidal cells. In Eq.@, the various quantities are given by 

L xy {l - (1 - x)[l - (1 - L xy )(l - y)}} 



Sl 



x + (1 - L xy )L xy (l - x) 2 (l - y) 



2 

—xy 

Fx 



2 ' 



3[x + (l-L xy )L xy (l-x) 2 (l-y)\ 
Np _ [L xy (l-x) + x}(l-x)(l-y) 

" 3[-x - (1 - L xy )L xy (l - x) 2 (l - y)} • 

Note that we have assumed x to be a real number, that is, x ~ o s jo2^ which will be 

justified below. After substituting e = e + a/i2irf into Eq.(|5|), we rewrite b xy after simple 

manipulations 

^ = (W P + _^_) + _| 7 _ (6 ) 

with s = (1 — ei/e2) _1 and t = (1 — ai/o^) -1 , where the dispersion magnitude and charac- 
teristic frequency admit respectively 

(t - si)(s - Sl) 

= 1 <r 2 s(t-si) 

2tt e 2 t(s - si)' 1 J 



III. NUMERICAL RESULT 



In Fig.l(a)-(c), the dependence of the pole, characteristic frequency and dispersion mag- 
nitude on the parameter x are investigated for different particle shape. For clarity, we set 
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z = 1/y from now on. It is evident that both s± and f c decreases monotonically as x in- 
creases; larger depolarization factor leads to smaller S\ and f c . However, for increasing x, 
Se increases concomitantly; also, Se increases while L z increases. 

In Fig.l(d)-(f), we investigate the dependence of the pole, characteristic frequency and 
dispersion magnitude on x for different shell thickness. For z — 1, namely without shell, 
si = 0.3(= L xy ) always. Increasing x leads to decreasing si or f c . It is evident that thicker 
shell yields larger pole or f c as x < 1, but smaller as x > 1. In the dispersion magnitude 
plot, Se is constant for the uncoated particles. In the case of coated particles, increasing x 
yields increasing Se; moreover, larger z, smaller Se. 

In Fig. 2, the dependence of material parameters (namely, s, t, e 2 , and <r 2 ) on f c and Se is 
discussed. f c is weakly dependent on s, but strongly dependent on t, e 2 and a 2 . Smaller \t\ 
(or e 2 ), larger f c . However, larger cr 2 , larger f c . On the other hand, Se is weakly dependent 
on s and t. Also, it is evident that e 2 and <r 2 play no role in Se, as Eq.(|7]) actually shows. 
Note all curves in Fig. 2(g) and (h) are respectively overlapped. In other words, this predicts 
that e 2 and <r 2 plays no role in the peak value of the angular velocity (see Fig.4). 

In Fig. 3, the dependence of L z , z and x on — Im[b xy ] are investigated, respectively. There 
is always a peak, the location of which is just the characteristic frequency f c . As L z (or 
x) increases, f c becomes red-shift, while the corresponding peak value increases as well. 
However, a larger z yields smaller peak value, while f c is changed weakly. 

In Fig.4, we discuss the dependence of s, t, e 2 and cr 2 on — Im[b xy }. The s effect is small 
so that it can be neglected. Smaller \t\ leads to larger f c ; at the same time, increasing peak 
value appears. Increasing e 2 yields decreasing f c , and the related peak value is unchanged. 
For increasing <r 2 , f c increases concomitantly, but the peak value remains unchanged as well. 

In Fig. 5, the experimental data is extracted from an experiment on human erythrocytes 
fl~4| . Obviously, there is a good agreement between theory and experiment. During this 
fitting, we model the cell as oblate spheroid [L z = 0.65). 
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IV. DISCUSSION AND CONCLUSION 



In the present work, we have discussed the effect of oblate spheroidal particle shape on 
the ER spectrum. In reality, there exist many cells in the form of oblate spheroid, such 
as human erythrocytes [[J. However, there exist only a few theories to discuss the ER 
spectrum of such cells (e.g., [13] and references therein). Our theory is advantageous in that 
the characteristic frequency at which the maximum rotational velocity occurs, is derived 
analytically, which simplifies the fitting of experimental data. 

We have assumed x, namely the ratio of the shell to host dielectric constant, to be a 
real and positive number. In fact, its imaginary part is indeed small. If we had retained the 
imaginary part of x, two peaks would appear, and the conductivity-dominated peak would 
have occurred at substantially lower frequency. Therefore, the neglect of the imaginary part 
of x is to drop the lower frequency peak. Moreover, according to our calculations, there 
are two (degenerate) sub-dominant poles associated with b xy in the spectral representation. 
Thus, in the present work, only one capacitance-dominated peak has been shown. 

In a recent work [f[D| , we developed simple equations to describe the ER of particles 
in a suspension from the spectral representation, and obtained a good fitting on a coated- 
bead ER assay by using the coated spherical model. In this work, an extension to oblate 
spheroidal case has been made. Our theory serves as a basis which describes the parameter 
dependence of the polarization and thereby enhances the applicability of various cell models 
for the analysis of the polarization mechanisms. In this connection, the shell-spheroidal cell 
model may readily be extended to multi-shell cell model. However, we have shown that the 



multi-shell nature of the cell may have a minor effect on the ER spectrum [17 



We have considered the isolated cell case, which is a valid assumption for low concen- 
tration of cells. However, for a higher concentration of cells, we should consider the mutual 
interaction between cells. When the volume fraction of the suspension becomes large, the 
particles may aggregate in the plane of the rotating applied field, and the mutual interactions 
between the suspended particles can modify the spin friction, which is a key to determine 
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the angular velocity of ER. For ER of two particles, we have successfully applied the spec- 
tral representation to deal with the dispersion frequency 0. However, the determination 
of the spin friction is still lacking. For two particles, the basic tool is the reflection method 
p~8[] , being analogous to the multiple image method in electrostatics, but being only valid 
for two particles. For more than two particles, we need a first-principles method, e.g., the 
Green's function (Oseen tensor) formulation. For a dilute suspension, however, one may 
adopt the effective medium theories |L9| to capture the effective viscosity of a suspension, 



hence modifying further the ER spectrum. 

One may argue that normal erythrocytes may deviate from an oblate spheroidal shape. 
For cells of non-conventional shapes, one may use a first-principles approach ||20|| , recently 
developed to deal with particles of rod shape. As shown in Ref. ||20|| , the derivation from 



spheroidal model can be small when the cells are rotating with their long axes along the 
applied field. 

In summary, we have presented a theoretical study of electrorotation based on the spec- 
tral representation theory. From the theoretical analysis, we find that cell shape as well as 
the coating can change the characteristic frequency. By adjusting the cell shape, dielectric 
properties and the thickness of the coating, it is possible to obtain good agreement between 
our theoretical predictions and the experimental data on human erythrocytes suspensions. 
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FIGURES 

FIG. 1. FIG.l. The dependences of pole, characteristic frequency and dispersion magnitude on 
x for s = 1.1, t = -0.005, e 2 = 80e , a 2 = 2.9 x 10 _5 5/m. (a)~(c) z = 2; (d)~(f) L z = 0.4. 

FIG. 2. FIG. 2. The dependence of characteristic frequency and dispersion magnitude on 
x for z = 6, L z = 0.4. (a) t = -0.005, e 2 = 80e , ct 2 = 2.9 x 10 _5 5/m; (b) s = 1.1, 
e 2 = 80e , a 2 = 2.9 x l^S/m; (c) i = -0.005, s = 1.1, ct 2 = 2.9 x 10~ 5 S/m; (d) s = 1.1, 
t = -0.005, e 2 = 80e ; (e) t = -0.005, e 2 = 80e , cr 2 = 2.9 x 10 _5 5/m; (f) s = 1.1, e 2 = 80e , 
a 2 = 2.9 x 10 _5 5/m; (g) 4 = -0.005, s = 1.1, ct 2 = 2.9 x 10 _5 5/m; (h) t = -0.005, s = 1.1, 
e 2 = 80e . 

FIG. 3. FIG. 3. The dependence of —Im[b xy ] on frequency for s = 1.1, t = -0.005. (a) 
e 2 = 80e , a 2 = 2.9 x 10 _5 5/m, z = 6, x = 2; (b) e 2 = 80e , ct 2 = 2.9 x 10~ 5 S/m, L z = 0.4, x = 2; 
( c ) e2 = 80e , ct 2 = 2.9 x 10 _5 5/m, z = 6, L z = 0.4. 

FIG. 4. FIG. 4. The dependence of —Im[b xy ] on frequency for x = 2, z = 6, L z = 0.4; (a) 
t = -0.005, e 2 = 80e , a 2 = 2.9 x 10 _5 5/m; (b) s = 1.1, e 2 = 80e , cr 2 = 2.9 x l^S/m; (c) 
t = -0.005, s = 1.1, a 2 = 2.9 x 10- 5 S/m; (d) t = -0.005, s = 1.1, e 2 = 80e . 

FIG. 5. FIG. 5. A fitting on human erythrocyte data. Experiment: a 2 = 12.5 x 10~ 3 . Theory: 
s = 100, t = 1.00783, e 2 = 80e , z = 3, x = 0.008, L z = 0.65, F/E 12 = 135. 
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